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Computer  Reduction  of  Meteorological  Rocketsonde  Data 

INTRODUCTION 


This  final  report  provides  complete  documentation  of  software,  and 
hardware  that  was  developed  to  permit  same  day  processing  of  meteoro¬ 
logical  rocketsonde  data  at  White  Sands  Missile  Range. 

Software  listings  of  all  programs  that  comprise  the  present  METROC 
program  are  included  along  with  detailed  flow  charts  for  those  programs 
that  were  specially  developed  to  handle  the  digitized  met  data.  Hardware 
documentation  includes  schematics  of  all  electrical  components,  overall 
diagrams  of  interconnection  wiring  and  the  microprocessor  computer  pro¬ 
gram  listing. 

As  an  indication  of  accuracy,  graphical  plots  of  significant  level 
temperature  profiles  as  determined  by  manual  methods  and  the  computer 
method  is  included. 

Since  the  report  that  covered  the  original  data  reduction  program'-’-' 
is  no  longer  available,  this  report  will  include  some  of  the  materia! 
from  that  report. 

Meteorological  rocket  sounding  systems  (metrockets)  have  been  devel¬ 
oped  '.o  obtain  upper  air  observations  in  the  lowest  100  km  of  the  atmo¬ 
sphere,  especially  that  portion  Inaccessible  to  routine  balloon  observa¬ 
tions,  i.e.,  above  30  km.  Some  of  these  systems  that  are  routinely  used 
include  a  payload  which  consists  of  an  atmospheric  sensor,  a  radio  trans¬ 
mitter  for  telemetering  the  measurements  to  ground  tracking,  receiving 
and  recording  equipment,  and  a  radar-reflective  retardation  device.  The 
atmospheric  sensor  normally  provides  temperature  measurements.  A  radar 

*A  General-Purpose  Meteorological  Rocket  Data  Reduction  Program  by  Mary 
Ann  Seagravcs,  ECOM  5  '2,  August  1972 
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track  of  the  payload  is  used  to  compute  wind  data  and  altitude  information 
to  correlate  with  the  temperature  measurements.  Another  t,pe  of  met  rocket 
payload  contains  only  a  radar  target,  i.e.,  chaff  or  inflatable  sphere, 
to  provide  wind  measurements.  Other  types  of  payloads  collect  specialized 
data  for  various  research  and  development  programs. 

The  purpose  of  the  computer  program  METROC  is  to  produce  routine 
wind  and  thermodynamic  data  from  the  raw  data  collected  from  metrocket 
firings  at  White  Sands  Missile  Range  (WSMR) ,  New  Mexico,  Utah  Launch  Com¬ 
plex  (ULC),  Green  River,  Utah,  and  other  sites. 

The  program  is  written  to  be  used  on  the  UNIVAC  1108  computer  at 
WSMR.  It  is  designed  so  that  the  complete  reduction  of  the  data  from 
each  firing  may  be  accomplished  with  one  pass  through  the  computer.  This 
increases  efficiency  in  operation  by  eliminating  delays  between  the  vari¬ 
ous  steps  of  the  data  processing.  With  this  approach,  however,  the  data 
may  not  be  monitored  during  intermediate  phases  of  the  reduction  process. 

The  data  input  to  METROC  may  be  both  winds  and  temperatures,  winds 
only,  or  temperatures  only.  A  digital  tape  of  the  radar  track  and  telem¬ 
etry  is  used  for  the  computation  of  wind  and  temperature  data  when  avail¬ 
able.  Manually  reduced  wind  and  temperature  data  may  be  input  using 
punch  cards. 

The  output  listings  are  in  metric  units  and  may  contain  any,  or  all, 
of  the  following: 

1.  Individual  wind  data  points 

2.  Summarized  data  for  significant  temperature  levels,  even  kilo¬ 
meter  levels,  5000  ft  levels  and  mandatory  levels 

3.  Listing  of  Meteorological  Rocket  Network  (MRN)  formatted  data 

4.  Ozonesonde  parameters 
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5.  Wind  data  in  English  units 

The  summarized  data  listings  may  include  layer  winds  averaged  uvc: 

2  kilometers,  corrected  temperatures,  pressure  data  computed  using  hydr  - 
static  equations  and  a  base-level  point  from  a  conjunctive  raoh,  atmo¬ 
spheric  density  and  speed  of  sound. 

Output  data  cards  are  punched  in  the  MRN  reporting  format  for  all 
runs  which  include  wind  data.  Additional  output  cards  are  punched  when 
the  English  units  output  option  is  selected. 

An  effort  has  been  made  to  conform  to  the  I R 1 G  standards  for  met- 
rocket  data  reduction  [1]  where  possible  throughout  the  data  reduction 
process.  Smithsonian  Meteorological  Tables  [2]  was  used  as  a  source  lor 
formulae  and  constants. 

1.  MATHEMATICAL  DEVELOPMENT 

A.  WIND  COMPUTATIONS  -  RADAR  DATA  TAPE  INPUT 

The  subroutine  READTP  obtains  the  raw  data  from  the  input 
tape  and  converts  it  to  standard  units  and  sampling  rates.  The  range 
is  converted  to  meters,  the  azimuth  and  elevation  to  radians,  and 
the  sampling  rate  to  10  samples  per  second.  The  sample  time  is 
referenced  to  time  of  launch.  The  raw  data  is  filtered  using  a 
technique  described  by  Avara  and  Miers  [J|.  A  117-point  symmetri¬ 
cal  filter  is  applied  to  data  so  that  5  consecutive  samples  of 
smoothed  data  are  available. 

Then,  for  each  of  these  5  points,  the  following  computations  are 

made : 

Altitude  corrected  for  earth's  curvature  Z 

2  _  R2  cos2  E _ 

c 

1.2742458  x  107  +  R  sin  E  +  7. 

A 


(1) 


where  R  =  range 

E  =  elevation  angle 

Z.  =  station  altitude. 

A 

Elevation  angle  corrected  for  earth's  curvature  F 

r  E  +  Z,  cos  E 

E  =  A 

c - 

1.40167038  x  107  sin  E. 

East-west  position,  X 

X  =  R  cos  E  sin  A 
c 

where  A  =  azimuth. 


(2) 


(3) 


North-south  position,  Y 

Y  =  R  cos  E^  cos  A.  (4) 

The  accelerations  are  then  used  to  correct  the  velocities  for 


sensor  motion  by  a  technique  developed  by  Eddy  [ 1 ] .  If  v  is  the  un¬ 
corrected  component  velocity,  then  v,  the  corrected  component  velocity, 


(5) 


where  a  =  component  acceleration 

v  =  vertical  velocity 
z  7 

a  =  vertical  acceleration 
z 


g  =  gravitational  acceleration. 

A  check  is  made  to  determine  if  apogee  has  been  reached.  When 
12  consecutive  points  have  been  encountered  for  which  the  altitudes 
are  decreasing,  it  is  assumed  that  apogee  has  been  reached  and  pay- 
load  expulsion  has  occurred. 

After  apogee  has  been  found,  additional  checks  are  made  to  deter 
mine  whether  the  data  is  valid.  If  the  computed  vertical  velocity 
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is  positive,  the  data  sample  is  rejected  as  invalid.  Other  (  hei  ks 

may  he  made,  hilt  there  is  an  option  to  hy-pass  these.  (See  sixth 

"plum  on  options  data  card).  They  are 

1  The  vert  11 al  velocity  must  he  greater  than  -250  n/s. 

2.  I  he  i  hangr  m  ti  ,al  vehnitv  from  the  previous  p>  •  1 1 1 1  di-rr 

•i  ’  exceed  10%. 

1  he  list  "t  ijues  t  i  onah  1  e  data  points  is  generated  that  i  ons  i  ■  '  : 
of  those  samp  It's  for  which 

1.  The  change  in  wind  direction  exceeds  30°  when  the  speed 
is  greater  than  20m/s,  or 

2.  There  is  a  100°/o  change  in  windspeed  when  the  wlndspecd  is 
greater  than  5  m/s. 

If  a  sample  is  found  to  be  invalid,  the  message  "DATA  EXCEEDS 
LIMITS  AT  TIMK  XXX"  is  printed  and  that  sample  is  not  used  in  con, lut¬ 
ing  layer  winds. 

When  a  sample  is  valid  (or  questionable),  it  is  printed  and  stored 
for  use  in  computing  layer  wind  data.  The  sample  output  rate  is  one  pet- 
second  . 

B.  WIND  COMPUTATIONS  -  CARD  INPUT 

The  wind  data  samples  input  on  punch  cards  arc  printed  and  stored 
for  subsequent  use  as  they  are  read.  A  message  is  printed  following 
each  point  for  which  the  time  or  altitude  is  not  sequential  or  for 
which  the  data  exceeds  certain  limits.  The  limits  Tor  valid  data  are 

0°  to  360°  for  wind  direction,  0  to  200  m/s  for  windspeed,  and  0  to 
250  m/s  for  fall  velocity. 

C.  COMPUTATION  OF  LAYER  WIND  DATA 

The  altitudes  of  the  points  for  which  layer  winds  are  to 
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be  computed  are  first  determined  in  the  main  program.  The  sub¬ 
routine  WNDAVG  is  called  to  compute  the  layer  wind  data. 

An  interval  over  which  averages  are  to  be  computed  is  deter¬ 
mined  for  each  layer  wind  data  point.  Then  each  of  the  individual 
data  points  is  read,  and  the  component  wind  velocities  are  accu¬ 
mulated  for  each  interval.  After  all  the  individual  data  points 
have  been  considered,  average  component  winds  and  fall  velocities 
are  calculated,  along  with  the  time  for  each  layer.  The  vector 
wind  is  determined  from  the  component  wind.  Any  layer  which  con¬ 
tains  fewer  than  n  points  is  bad-flagged,  n  being  6  when  the  wind 
data  is  computed  from  a  radar  tape  and  n=l  when  the  wind  data  is 
input  on  cards. 

D.  BASE-LKVEL  POINT 

A  base-level  point  is  a  data  point  obtained  from  a  conjunc¬ 
tive  rawinsonde  observation  that  was  taken  close  in  time  and 
space  to  the  rocketsonde  observation.  The  data  provided  is  at¬ 
mospheric  pressure  and  temperature  at  a  specified  altitude.  By 
using  this  and  rocket  temperatures,  thermodynamic  data  is  computed 
with  the  hydrostatic  equations.  To  be  valid,  the  temperatures 
measured  by  the  rawisonde  and  rocketsonde  should  agree  within  2.5° 
at  the  base-level  point.  The  program  will  reject  all 
base-level  points  that  do  not  have  this  agreement,  unless  a  special 
option  is  selected. 

Any  run  which  contains  temperature  data  may  contain  zero  to  5 
hase-level  points.  If  more  than  one  base-level  point  is  input,  the 
one  which  is  closest  to  25  kilometers  and  which  meets  the  tempera¬ 
ture  agreement  criterion  is  used  in  the  thermodynamic  computations. 
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Thermodynamic  computations  are  not  made  if  no  base-level  point 
is  input  or  if  none  meets  the  temperature  agreement  criterion  and 
the  option  overriding  this  criterion  is  not  selected.  These  com¬ 
putations  are  not  done  if  there  is  no  valid  base-level  point  avail¬ 
able  which  is  within  2  kilometers  of  the  lowest  temperature  data 
point. 

If  the  override  option  is  selected  and  none  of  the  input  base- 
level  points  meets  the  temperature  agreement  criterion,  that  point 
which  has  the  best  temperature  agreement  is  used. 

The  base-level  point  which  is  ,.o  be  used  in  thermodynamic 
computations  is  stored,  and  the  rocket  temperature  in  degrees 
Kelvin  and  the  geopotential  altitude  are  found: 

Tfc  *  T  273.16  (6) 

where  =  temperature  °K 
T  =  temperature  *C 

and 

zp  =  g  *  Re  x  Z  (7) 

9.8  (Re  +  Z) 

where  Z ^  a  geopotential  altitude 

g  =  gravitational  acceleration 

Re  a  local  radius  of  earth 
Z  a  geometric  altitude. 

E.  TEMPERATURE  AND  THERMODYNAMIC  DATA 

At  each  step  in  the  program  where  temperature  and  thermo¬ 
dynamic  data  are  to  be  processed,  the  altitudes  of  the  points  to 
be  output  are  determined  first.  The  altitudes  for  the  significant 
levels  are  interpolated  from  the  individual  wind  data  points; 
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when  wind  data  is  not  being  processed,  the  altitudes  are  read  from 
the  input  data  cards.  Significant  levels  are  not  processed  lor 
those  points  at  which  the  altitudes  are  not  available  from  the 
wind  data. 

Constant-altitude  layers,  i.e.,  even  kilometers,  0.5  kilo¬ 
meters,  500  feet,  or  1000  feet,  are  determined  by  considering  the 
highest  and  lowest  available  data,  either  wind  data  or  signifi¬ 
cant  temperature  levels.  A  table  of  constant-altitude  layers  is 
produced  for  points  within  this  span. 

A  table  of  standard  constant  pressures  is  contained  within 
the  program.  It  is  determined  which  of  these  pressures  fall  with¬ 
in  the  span  of  available  data,  and  then  altitudes  are  exponential¬ 
ly  interpolated  for  the  constant  pressure  (mandatory)  levels  from 
the  altitude-vs . -pressure  data  for  the  significant  levels. 

When  wind  data  processing  is  to  be  included,  the  layer  winds 
and  fall  velocity  are  computed  after  determining  the  layer  altitudes. 

Next,  except  for  significant  levels,  temperatures  for  each  layer 
are  linearly  interpolated.  Subroutine  THERMO  is  then  called  to  do 
the  thermodynamic  computations. 

First,  the  absolute  temperatures  are  computed  as  in  Equation 

(6). 

Then,  if  a  temperature  correction  is  to  be  computed,  a 
correction  is  derived  based  on  a  method  developed  by  Krumins 
and  Lyons  and  is  the  IRIG  temperature  correction  [4|: 

AT(k)  =  T(k)  -  A(z)  x  (Vz(k))2  +  B-^z)- x  (T(-k ~-T(k~1) 

t (k+ 1 )  -  t(k-l) 

-C(z)  +  D(z)  x  T(k)4  (8) 
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where  AT(k)  =  temperature  correction  for  kth  point 
T(k)  =  absolute  temperature  at  kth  point 
t(k)  =  time  of  kth  point 
Vz(k)  =  fall  velocity  at  kth  point 
and  A(z),  B(z),  C(z)  and  D(z)  are  coefficients  which  vary  with 
sensor  type  and  altitude  and  C(z)  is  also  a  function  of  time  of  day 
(night  or  day). 

NOTE:  When  wind  data  is  not  processed  in  a  run,  fall  velocity 
data  is  not  available.  Thus,  no  provision  has  been  made  to  compute 
temperature  corrections.  However,  temperature  corrections  which  have 
been  manually  computed  may  be  input  on  data  cards  and  used  in  the  final 
computat ions . 

The  temperature  correction  is  added  to  the  uncorrected  temperatures, 
and  the  geopotential  altitude  is  computed  for  each  point  as  in  Equation 
(7). 


Atmospheric  pressure  is  computed  next,  except  when  processing  con¬ 
stant-press  levels.  Any  time  there  is  a  layer  of  missing  temperature 
data  greater  than  3  kilometers,  pressure  and  related  computations  are 
not  continued  beyond  that  point.  Pressure  is  computed  both  upward  and 
downward  from  the  base- level  point: 

7  Clrl 

(9) 


P(k)  =  P  (k-1)  x  exp  ^p  (k-1)  ^p^ 


14.63725  x  (T(k)  -  T(k-l)) 

where  P(k)  =  pressure  (mb) 

Zp(k)  =  geopotential  altitude  (m) 

T(k)  =  corrected  absolute  temperature 
k-1  is  the  previous  point  when  k> I ,  and  k-1  is  the  base-! eve 
point  when  k=l. 
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Density  is  computed: 

D  =  P  x  348.38/T  (10) 

D  =  D  x  1.9403806  x  10'6  (II) 

where  D  =  density  (gm/cu  m) 

P  =  pressure  (mb) 

T  =  corrected  absolute  temperature 
Dj  =  density  (slugs/cu  ft). 

Speed  of  sound  is  calculated: 

S  =  20.0544  x  VT  (12) 

where  S  =  speed  of  sound  (m/s) 

T  =  absolute  temperature 
I.  SYSTEM  FORMULATION 

In  order  to  process  high  altitude  meteorological  data  from  a 
rocket  born  datasonde,  the  radar  positional  data  and  the  significant 
level  temperature  data  from  the  datasonde  must  be  combined  into  a 
single  prof i le . 

The  system  described  in  this  report  provides  a  means  of  com¬ 
bining  both  the  radar  data  and  telemetry  data  into  a  common  data 
bank  for  final  processing  by  the  UNIVAC  1108  computer  into  a  meteo¬ 
rological  profile. 

A.  OUTLINE  OF  TOTAL  SYSTEM 

The  total  system  consists  of  l)-a  microprocessor  system  that 
will  digitize  the  temperature  telemetry,  form  the  appropriate  mes¬ 
sage  format  with  this  data,  and  transmit  this  data  to  the  control 
record  facility  -  this  combines  the  temperature  data  but  not  signif¬ 
icant  level  data  with  the  radar  data;  2)-software  or  that  will 
process  the  digitized  telemetry  data  to  form  a  significant  level 
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temperature  profile  that  is  consistant  with  the  levels  selected  by 
hand.  (4) 

Once  the  data  is  in  significant  level  form,  the  final  existing 
computer  program  (ROCKET)  will  generate  the  final  met  profile. 

B.  MICROPROCESSOR  SYSTEM 

The  microprocessor  chosen  was  the  INTEL  8080  configured  as  an 
SBC  80  single  hoard  computer.  [  5  J  This  provides  the  user  with  1 02* 
bytes  of  RAM,  1024  bytes  of  ROM  or  EPROM,  two  8255  parallel  1/0,  ports, 
one  8251  USART  serial  1/0  port,  EIA  RS  232  logical  level  converters, 
card  cage  and  power  supply.  This  particular  system  was  selected 
because  it  has  the  necessary  hardware  to  do  the  job  and  was  available. 

C.  ADDITIONAL  SYSTEM  REQUIREMENTS 

In  order  to  provide  additional  future  capabilities,  provisions 
were  made  to  read  and  incorporate  into  the  message,  time  data  from  an 
IRIG  B  time  code  reader,  and  data  from  four  event-marker  flip-flops. 

These  hardware  items  together  with  the  analog-to-digi tal  con¬ 
verter,  and  tri-state  multiplexing  chips  were  mounted  on  a  general 
purpose  prototype  board.  This  board  was  inserted  in  the  card  cage 
to  form  the  total  system.  A  block  diagram  of  this  system  is  shown 
in  figure  1  . 

D.  MESSAGE  FORMAT 

The  message  format  for  transmission  to  the  Central  Record 
Facility  (CRF)  is  composed  of  two  120  hit  messages.  The  only 
restrictions  in  these  messages  is  that  the  first  15  hits  are  Fixed 
by  the  sync  and  station  identification  numbers.  The  message'  format 
that  was  used  is  shown  in  table  1.  In  looking  at  table  1  it  should 
be  noted  that  the  time,  seconds,  is  contained  in  two  seperale  bit 
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Figure  1.  System  block  diagram 
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blocks.  This  problem  occurred  because  the  original  time  code 
reader,  for  which  the  system  was  designed,  had  time  broken  down 
into  BCD  10  h,  hours,  10  min,  min,  10  s,  and  sec.  When 
the  modems  were  installed,  a  new  time  code  reader  was  also  installed 
that  had  time  in  binary  seconds  and  milliseconds  as  indicated. 

The  new  data  was  brought  into  the  system  by  changing  only  the  wiring 
to  the  connector.  The  time  data,  met  data,  and  event-marker  data  is 
combined  properly  by  the  software  program  that  generates  the  signifi¬ 
cant  level  data. 

K.  8080  SOFTWARE 

Figure  2  is  a  flow-chart  of  the  software  program  that  performs 
the  data  collection,  message  formation  and  transmission.  The  last 
block  in  figure  2  was  necessary  because  the  data  when  received  by 
the  CRF  is  complimented  and  the  1st  bit  received  is  treated  as  the 
most  significant  bit  of  the  message.  In  contrast,  the  modem  Interface 
(USART)  sends  the  least  significant  bit  first.  The  software  program 
for  the  ROM  control  corresponding  to  figure  2  is  shown  in  table  2. 

III.  MICROPROCESSOR  SYSTEM  HARDWARE  DIAGRAMS 

The  system  is  composed  of  six  subsystems.  1)  -  SBC  80  micro¬ 
processors,  2)  -  Interface  card,  3)  -  I R I G  B  time  code  reader-gen¬ 
erator,  A)  -  Modem,  5)  -  Card  cage,  6)  -  Power  Supply. 

Items  1,  3,  4,  5  and  6  are  standard  commercial  units  and  will 
not  be  described  in  detail. 

A.  INTERCONNECTION 

Figure  3  illustrates  the  various  connection  cables  that  are 
used  to  interconnect  the  various  subsystems.  The  only  connection 
not  indicated  is  the  backplane  connection  -  this  is  made  when  the 
two  cards  are  inserted  into  the  card  cage. 
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Flow  chart  for  8080  software. 
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Figure  2 
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TABLE  Z  (cont) 
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TABLE  2  (cont) 


PROGRAM  IS  223  BYTES  LONG  (1)  WITH  0  ERRORS  DETECTED. 

L  SYMBOL  LISTING  ( Y=YES , N=NO) ? 

CLRW=002B  DOIT=00C4  G0=0020  HUSH=000D  ID=0053  MODE=0001  NEWW=0030 
NOGO=001 1  PSW=0006  RSET=0000  SP=000f-  SWAP=00CE  SYN1=0049  SYNC=0005 
WD10=009D  WD1 1 =00A7  WD16=00AC  WDl7=00B4  WD18=00BC  WORD=3COO  WRD1=0044 
WRD2=004E  WRD3=0058  WRD4=0061  WRD5=006B  WRD6=0075  WRD7=007F  WRD8=0089 
WRD9=009 3  XMIT=00C8  ZERO=005D 


— 


Table  3  shows  the  pin  connections  and  signal  identification  lor 
all  connectors.  TB  refers  to  Terminal  Block. 

B.  INTERFACE  CARD 

Figure  4  shows  the  component  locations  on  the  interface  card. 
Each  integrated  circuit  (IC)  is  indicated  by  a  number  (1  thru  17) 
and  IC  type  (SN  74367,  etc.). 

Figure  5  illustrates  the  event  marker  logic,  Go  and  Reset  logic 
and  the  BO  thru  B3  data  flow  from  the  Analog  to  Digital  Converter 
(A/D).  The  8  bit  word  composed  of  B3,  B2,  B1 ,  BO  and  event  marker 
bits  1  thru  4  are  gated  into  the  microprocessor  when  pin  9,  IC  13  is 
active.  The  numbers  inside  the  gates  indicate  the  1C  number.  The 
event  marker  outputs  are  set  by  a  +5  volt  pulse. 

Figures  6  thru  8  show  the  various  bits  that  compose  the  words 
brought  into  the  microcomputer.  Figure  9  illustrates  the  A/D 
pin  connections  and  the  multiplexer  circuit  for  selecting  the  various 
words . 

An  ADC  (Analog-to-Digital  Converter)  manufactured  by 
Burr-Brown  Research  Corporation  was  used  to  digitize  the 

met  data.  The  met  data  signal  was  available  from  an  on-site  fre¬ 
quency  to  voltage  converter. 

IV.  SOFTWARE 

A.  ORIGINAL  METHODS 

In  order  to  process  high  altitude  meteorological  data  from  a 
rocket  born  datasonde,  the  positional  data  from  the  tracking  radar 
and  the  significant  level  temperature  data  from  the  datasonde  must 
be  combined  into  a  single  profile. 
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Figure  4.  Component  location  on  Interface  card. 
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Figure  5.  Event  Markers,  Go,  Reset  Diagram. 
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Originally  the  telemetry  data  was  processed  by  hand  to  form 
significant  level  temperature  data.  This  data  is  a  table  oi  time 
after  launch,  and  corresponding  temperatures,  which  were  punched  on 
cards  for  the  later,  combining  with  the  radar  positional  data. 

The  radar  positional  data  of  azimuth,  elevation,  and  slant  range 
of  the  datasonde  are  transmitted  via  a  syncronous  modem  to  a  central 
record  facility  (CRF)  and  recorded  on  a  uni  log  magnetic  tape.  Later 
the  ROCKET  program  calculates  a  wind  profile  from  the  radar  data 
and  combines  the  temperature  data  with  it  from  the  significant  tem¬ 
perature  cards  -  to  produce  the  final  met  profile. 

B.  REAL-TIME  METHOD 

The  real-time  system  records  simultaneously  the  digitized  met 
data  from  the  microprocessor  system  and  the  positional  radar 
data  on  uni  log  magnetic  tape. 

When  the  unilog  tape  has  been  transferred  to  the  1108  system 
B  computer  system,  the  modified  ROCKET  program  then  processes  both 
sets  of  data  to  form  significant  level  temperature  data,  radar 
wind  data,  and  ultimately  a  final  met  profile. 

V.  DESCRIPTION  OF  NEW  SOFTWARE 

In  order  to  process  the  telemetry  data  to  get  significant 
level  temperature  data,  6  new  programs  were  developed.  A  brief 
description  cf  these  programs  and  corresponding  flow  charts  follow: 
Subroutine  SIGLEV  I1SIG,  TIME,  VALUE,  REFT) 

As  the  telemetry  is  read  from  the  UN1L0G  tape  by  RTDATA,  SIGLEV 
is  called  to  build  an  array  of  reference  ordinates  and  significant 
temperature  ordinates,  which  have  passed  through  the  filter,  along  with 
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their  corresponding  times.  These  ar-  passed  to  the  subroutine  by  the 
arguments  TIME  and  VALUE.  The  significant  temperature  ordinates  are 
selected  by  a  call  to  subroutine  SIGPT. 

When  the  entire  tape  is  read,  the  significant  temperature  ordi¬ 
nates  are  converted  into  ordinate-to-reference  ratios.  Once  an 
ordinate-to-reference  ratio  is  interpolated,  subroutine  LoKTEM  con¬ 
verts  the  ratio  to  an  actual  temperature. 

Finally,  after  the  conversion  of  the  data  to  a  temperature  array, 
the  significant  levels  of  this  array  are  found  within  a  tolerance  ot 
2°C  by  subroutine  SIGPT.  The  final  significant  temperatures  art'  then 
printed  along  with  their  corresponding  flight  times,  ratios,  and  ordi¬ 
nate  values. 

The  other  two  calling  arguments  do  the  following:  11S1G  flags 
the  subroutine  for  initialization  and  final  runs,  and  HEFT  is  the 
distinguishing  level  for  testing  whether  VALUE  is  a  reference  or  a 
temperature  ordinate. 

S ubrout ine  LOKTEM  (1CAL,  CRATIO) 

The  subroutine  LOKTEM  h3s  two  functions  which  are  distinguished 
by  the  argument  I CAL. 

The  first  is  performed  when  RTDATA  calls  LOKTEM  to  read  temper¬ 
ature  calibration  data  from  cards  and  to  set  up  an  interpolation  scheme 
for  temperature  calculation.  Once  this  is  done,  a  check  of  the  inter¬ 
polation  scheme  is  made.  If  the  calculated  values  differ  from  the  cali¬ 
bration  values  by  more  than  1°C  the  baseline  calibration  is  rejected 
and  a  message  written.  The  calibration  ratios,  values,  and  test  are  all 
printed  on  the  output  as  well. 


-29- 


The  second  function  is  to  perform  normal  temperature  calculations. 

A  temperature  ordinate  to  reference  ordinate  ratio  (CRAT10)  is  passed 
to  the  subroutine  by  SIGLEV  as  it  builds  the  arrays.  From  this  ratio 
the  subroutine  interpolates  a  temperature  value  from  the  array  of  stored 
calibration  values  and  returns  it  via  CRAT10. 

Subroutine  RTDATA 

RTDATA  is  a  general  purpose  tape-read  routine  used  when  processing 
data  from  a  UN1L0G  tape.  Subroutine  READTP  and  TAPE1  use  RTDATA  to 
read  FPS-16  radar  data,  as  well  as  MET  data  for  processing  temperatures 
from  telemetry.  When  RTDATA  is  called,  the  data  in  all  active  sub¬ 
channels  for  one  time  is  stored  in  the  COMMON  area  FPS.  Since  the  data 
is  recorded  on  the  tape  in  blocks  of  20  samples  per  record,  a  new  record 
is  read  once  every  20  times  the  subroutine  is  called. 

If  temperature  telemetry  is  being  processed,  RTDATA  initially  calls 
LOKTEM  to  perform  the  baseline  temperature  calibration.  During  subse¬ 
quent  calls,  the  data  is  filtered  by  calling  subroutine  FILTER  and 
the  calculation  of  temperatures  is  overseen  by  calls  to  SIGLEV. 

When  an  EOF  indication  is  read,  the  final  call  to  SIGLEV  is  made  to 
finish  the  calculations  and  print  the  significant  temperature  levels. 

When  no  temperature  telemetry  is  processed  from  the  tape,  radar- 
data  is  read  as  needed  until  an  EOF  occurs. 

Subroutine  SIGPT 

This  subroutine  determines  which  temperature  out  of  a  string  are 
significant  -  namely  those  points  which  when  connected  by  straight 
lines  forms  a  temperature  profile  that  is  within  a  given  tolerance  of 
the  true  profile. 
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Each  temperature  ordinate  is  taken  as  it  comes  and  an  array  of 
significant  temperatures  is  formed.  After  processing  the  ordinate 
values,  the  significant  points  are  converted  to  degrees  and  passed 
through  the  subroutine  a  second  time  with  the  tolerance  set  at  two 
degrees . 

S u i >r°ut_ine  FI  [.TER 

This  subroutine  creates  a  fi 1 tered-ordinate  sum  of  six  values. 
This  is  done  by  first  selecting  every  eighth  point,  and  then  comparing 
each  point  to  adjacent  points  for  correlation.  When  an  ordinate  value 
survives  the  filter,  it  is  added  to  the  previous  sum  until  a  string  of 
six  is  completed.  The  sum  is  passed  to  the  calling  program  and  a  flag 
is  set  to  indicate  that  a  filtered  ordinate  value  is  available. 

Sub routine  R0C0B1 

The  subroutine  R0C0B1  is  called  by  the  main  program  to  perform 
the  task  of  formulating  and  printing  the  ROCOB  message. 

R0C0B1  is  sent  the  complete  array  of  1-KM  level  wind  and  thermo¬ 
dynamic  data  in  the  COMMON  area  ROCOB. 

VI.  DESCRIPTION  OF  REMAINING  SUBROUTINES 
Sub roii t i ne  MRN  ( NOM ,  ID,  IV  IC; 

The  subroutine  MRN  arranges  the  wind  and  thermodynamic  data 
found  in  the  COMMON  areas  WIND  and  TMP  into  the  Meteorological 
Rocket  Network  reporting  format  [A],  prints  it,  and  punches  it 
on  data  cards.  Since  leading  zeros  are  to  be  preserved  in  this 
format,  all  interger  fields  are  broken  into  fields  of  one  di<git 
so  that  zeros  will  appear  when  printed  and  punched  via  Fortran 
statements . 
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The  variable  NOM  is  the  number  of  points  to  be*  processed. 

Tht*  variable  ID  equals  30  when  processing  significant  or  kilo¬ 
meter  levels,  40  for  constant-pressure  levels.  The  variable  1 V. 
is  an  indicator  of  wind  input  type,  while  IC  indicates  uncorrected 
or  corrected  winds . 

The  COMMON'  area  MRN  is  used  to  store  the  data  as  it  is  broken 
into  single-digit  fields.  The  COMMON  area  SIGN  is  made  up  of  alpha¬ 
numeric  fields  containing  a  plus  sign,  a  minus  sign,  a  blank,  and 
a  n i ne  . 

Subroutine  OZONE  (Nl,  N2) 

The  subroutine  OZONE  computes  and  prints  special  thermody¬ 
namic  parameters  which  are  of  significance  in  reference  to  met- 
rockets  with  ozonesonde  payloads.  These  parameters  are  d  (In  (T)) 
/dZ  and  d(ln(P))/dZ,  where  T  is  temperature,  P  is  pressure,  and  Z 
is  altitude. 

The  data  used  in  the  computations  is  found  in  the  COMMON  areas 
WIND  and  TMP.  The  variables  Nl  and  N2  are  the  indices  of  the  first 
and  last  points  to  be  used. 

Si u b routine  QUESS 

The  subroutine  QUESS  prints  the  time  and  altitude  of  as  many 
as  100  wind  data  samples  of  questionable  validity.  A  point  is  con¬ 
sidered  questionable  if  either  there  is  a  direction  shift  of  at 
least  308  when  the  windspeed  Is  greater  than  20  m/s,  or  there  is 
a  100  percent  change  In  speed  with  a  windspeed  greater  than  5  m/s. 

The  times  and  altitudes  of  the  questionable  points  are  stored 
by  the  main  program  In  the  COMMON  area  QQ,  along  with  the  number  of 
questionable  points. 


-32- 


S u b r 0 u t  in e  TPREAD  (10F,  I ND ,  NEN'D) 


The  subroutine  TPREAD  is  used  to  obtain  raw  data  from  digital 
tape  and  convert  the  times,  units,  and  sampling  rates  to  standard 
values.  For  each  sample,  the  data  is  stored  in  the  COMMON  area 
IRAK.  Time  of  the  sample,  TM,  is  in  seconds  referenced  to  lift-off 
time.  Range,  RG,  is  in  meters  while  azimuth,  AZ ,  and  elevation, 

EL,  are  in  radians.  The  standard  sampling  rate  is  10  samples  per 
second.  An  end  of  file  indicator  is  found  in  the  parameter  10F. 

The  variable  IND  specifies  the  type  of  data  tape  being  processed. 
NEND  contains  the  number  of  reels  of  input  tape. 

When  the  input  radar  data  is  contained  on  a  DR  format  tape, 
the  time  and  sampling  rate  are  not  altered.  The  range  is  converted 
from  yards  to  meters,  and  the  azimuth  and  elevation  are  converted 
from  mils  to  radians. 

If  the  input  data  is  on  a  UNILOG  tape,  the  sampling  rate  is 
reduced  from  20  samples  per  second  to  10  samples  per  second  by 
averaging  each  2  consecutive  points.  The  time  is  referenced  to  the 
time  of  launch  by  subtracting  the  lift-off  time,  found  in  the  COMMON 
area  LIFTTM,  which  also  contains  input  subchannel  numbers  and  radar 
numbers.  The  range  is  converted  from  yards  to  meters,  and  the  an¬ 
gles  are  converted  to  radians.  If  two  radars  are  to  be  processed  in 
one  run,  the  data  for  the  second  radar  is  converted  and  written  on 
Fortran  logical  unit  3,  while  the  data  for  the  first  radar  Is  being 
reduced.  This  data  is  used  directly,  without  further  conversion, 
when  processing  the  second  radar. 

The  data  is  read  from  the  UNILOG  tape  by  the  subroutine  RTDATA 
and  is  stored  in  the  common  area  FPS.  The  common  area  EOFF  contains 
an  end  of  file  indicator. 
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Subroutine  THERMO  (TIEP,  ALT  IT,  TIET  Nl,  N2,  NON,  IN.NCOR,  GG,  RE) 
The  subroutine  THERMO  computes  corrected  temperature,  pressure, 
density  and  speed  of  sound  for  the  data  points  found  in  the  COMMON 
area  WIND  and  stores  them  in  the  COMMON  area  TMP.  The  uncorrected 
temperature  corresponding  to  each  point  in  the  array  ARR  is  stored 
in  the  array  ARC  in  the  COMMON  area  TMP. 

First,  t he  temperature  correction  is  computed  by  using 

Equation  (8).  The  geopotential  altitude  is  computed  thus: 

z  _  g  x  Re  x  Z 

P  - - - 

9.8  x  (Re  +  Z) 

where  Z^  =  geopotential  altitude 

g  =  acceleration  due  to  gravity 

R  =  radius  of  the  earth 
e 

Z  =  geometric  altitude. 

Except  when  processing  constant-pressure-level  data,  the  next 
step  is  to  compute  corresponding  pressure.  For  the  kth  point, 

P(k)  =  P(k-l)  x  exp  (Z  (k-l)-Z  (k)/(14.63725x(T(k)+T(k-l)))| 
where  P  =  pressure  (mb) 

Zp  =  geopotential  altitude  (m) 

T  =  temperature  (°K) 

and  k-1  is  the  previous  point  computed  if  k>l,  the  raob  base- 
level  point  is  k=l. 

Then  density  is  computed: 

D  =  P  x  348.38/T 
where  D  =  density  (gm/cu  m) 

P  =  pressure  (mb) 

T  =  temperature  (°K) . 
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Also, 


Dj  =  D  x  1.9403806  x  10"6 
where  =  density  (slugs/cu  ft). 

Speed  of  sound  is  calculated: 

S  =  20.0544  X  VT 
where  S  =  speed  of  sound  (m/s) 

T  =  temperature  (°K) . 

All  missing  data  is  bad-flagged. 

The  variables  TIEP,  ALTIT,  and  TIET  are  the  RAOR  base-level 
pressure,  geopotential  altitude  and  rocket  absolute  temperature. 

The  variables  N1  and  N2  are  the  indices  of  the  first  and  last  points 
for  which  there  is  temperature  data,  while  NON  is  the  number  of 
points  in  the  array  ARR.  The  variable  IN'  indicates  the  type  of  data 
levels.  The  temperature  correction  to  be  used  is  specified  by  the 
variable  NCOR.  GG  is  the  local  acceleration  due  to  gravity  and  RE 
is  the  local  radius  of  the  earth. 

Subroutine  WNDAVG  (DIV,  NOM ,  KK,  JL) 

The  subroutine  WNDAVG  computes  averaged  wind  data  for  an 
interval  of  size  DIV.  Thc__wind  data  samples  which  are  to  he 
used  in  averaging  were  written  on  Fortran  logical  unit  13  as 
they  were  computed  by  the  main  program.  As  the  individual  sam¬ 
ples  are  read,  sums  of  the  x,  y,  and  z  component  velocities  are 
computed.  After  all  the  data  has  been  read,  component  velocities 
are  computed.  Vector  wind  and  wind  shear  are  then  calculated 
from  the  component  velocities.  The  time  is  obtained  by  interpola¬ 
tion.  Any  layer  for  which  fewer  than  N  points  wore  found  is  bad- 
flagged,  where  N  is  6  when  the  wind  is  computed  from  radar  input 


A 
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tape,  and  N  equals  1  when  the  wind  is  input  on  data  cards. 

All  computed  data  is  stored  in  the  array  ARR  in  the  COMMON 
area  WIND.  The  variable  NOM  is  the  number  of  points  for  which 
wind  averages  are  to  be  computed.  The  variable  KK  indicates 
whether  the  altitudes  to  be  used  are  in  meters  or  feet.  The 
variable  JL  is  used  to  indicate  the  wind  data  input  type,  i.e., 
cards  or  digital  tape. 

VII.  GLOSSARY  OF  MNEMONICS 


Main  Program 

A 

AA 

AB 

ACC  Cl, 10) 

ALT  (200) 

ALT  IT 
ALT I TD 
APO  (12) 

ARC  (7,200) 
elt  1 
2 
3 
A 

5 

6 
7 

ARR  (10,200) 
elt  1 
2 
3 
A 

5 

6 

7 

8 
9 

10 

ATE 

AZ 

BLANK 

BLANN 


wind  direction  in  radians 

constant  used  in  computing  constant- 

pressure  altitudes 

constant  used  in  interpolating 

array  used  in  storing  acceleration  prior 

to  use  of  10-point  filter 

array  containing  altitudes 

geopotenti.il  altitude  of  base-level  point 

geometric  altitude  of  base-level  point 

array  of  altitudes  used  in  searching  for 

apogee 

array  used  to  store  tbermodynami c  data 

temperature  °C 

temperature  °K 

temperature  correction 

pressure  mb) 

density  (gm/m3) 

density  (slug/cu  ft) 

speed  of  sound  (m/s) 

array  used  to  store  layer  wind  data 

time  (sec) 

altitude  (meters) 

altitude  (feet) 

N-S  wind  component 
E-W  wind  component 
windspeed 
wind  direction 
fall  velocity 
wind  shear 
time  (min  sec) 

layer  thickness  used  in  computing  altitude 

array 

azimuth 

5-character  field  of  blanks 
1-characler  field  containing  a  blank 
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c 

CONF 


CONM 

CONI 

CONZ 

DF.N 

DY 

EL 

FLAG  (200) 

GG 

IDD 

IND 


INK 
I  OF 

I  OF  (11) 
olt  1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

IRAD  (2) 
l  ROAR 
ISC  (2) 
1TE 

K 

KEX 

KK 

KLM 

KYM 

LAST 

LF1RST 

LI. 

MINUS 


constant  equal  to  .0174533 

factor  used  to  convert  input  altitudes 

to  feet 

factor  used  to  convert  input  altitudes 
to  meters 

1.94254  converts  m/s  to  knots 
-1.68895  converts  knots  to  ft /sec- 
constant  equa I  to  change  in  time 
day  of  month 
elevation 

array  used  to  indicate  missing  tempera¬ 
ture  data 

acceleration  due  to  gravity 
data  card  type  indicator 
indicates  type  of  data  being  processed 
=  -1  when  processing  one  radar  from 
UN  I  LOG  tape 

=  0  when  processing  1st  of  2  radars 
(UNI  LOG  tape) 

=  1  when  processing  2nd  of  2  radars 
(UNI  LOG  tape) 

=2  when  processing  DR  format  tape 

indicator  =  0  except  when  missing  data 

layer  greater  than  5  km 

end  of  file  indicator 

input  options 

wind  input  type 

temperature  input  type 

temperature  correction  type 

English  units  output  option 

ozone  parameters  output  option 

constraints  on  wind  data  option 

corrected  winds  input  indicator 

use  base-level-point  option 

continue  temperature  data  indicator 

continue  thermodynamic  data  option 

suppress  printing  1/sec  wind  data  option 

radar  number 

radar  being  processed 

input  subchannel 

layer  thickness  used  in  computing  alti¬ 
tude  array 

temperature  correction  type  indicator 
number  of  questionable  data  points 
index  used  in  searching  for  apogee 
=  0  until  printing  wind  shear  data 
=  0  until  first  wind  data  card  has 
been  checked 

altitude  of  last  point  in  array 
altitude  of  first  point  in  array 
index  used  in  Eddy  correction 
1-character  field  containing  a  minus 
sign 
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MK  (11) 

elL  1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

MON 

MONTH  (24) 
MS1G 
MUN  (2) 

Ml  (5) 

M2  (3) 

M3  (3) 

M4  (3) 

M5  (3) 

M6  (3) 

M7  (3) 

M8  (2) 

Mil  (3) 

NEND 

NINE 

NMP 

NOM 

NOMM 

NOMP 

NOM  I 

NZ 

NI 

N2 

PLUS 

QUS  (2,100) 

RAW  (4,121) 
elt  1 
2 

3 

4 

RDNUM 

RE 

RG 

ROCKT  (3) 
SALT 


array  containing  data  to  be  output  in 

MRN  format 

altitude 

wind  direction 

windspeed 

N-S  wind  component 
E-W  wind  component 
fall  velocity 
temperature 

temperature  correction 

power  of  10  multiplied  to  obtain  pressure 
power  of  10  multiplied  to  obtain  density 
speed  of  sound 
month  (integer) 

array  containing  names  of  months 
number  of  base-level  points  input 
two  digits  representing  month 
five  digits  representing  altitude 
three  digits  representing  wind  direction 
three  digits  representing  windspeed 
three  digits  representing  N-S  wind  com¬ 
ponent 

three  digits  representing  E-W  wind  com¬ 
ponent 

three  digits  representing  fall  velocity 
three  digits  represent ing  temperature 
two  digits  representing  temperature 
correction 

three  digits  representing  speed  of  sound 
number  of  reels  of  input  tape 
1-character  field  containing  a  nine 
number  of  significant  temperature  points 
input 

number  of  levels  in  array 
number  of  levels  in  array 
number  of  levels  in  array 
number  of  levels  in  array 
number  of  consecutive  bad  wind  data  points 
index  of  first  level  with  temperature  data 
index  of  last  level  with  temperature  data 
1-character  field  containing  a  plus  sign 
array  containing  time  and  altitude  ol  ques¬ 
tionable  wind  data  points 
raw  radar  data  array 
time 
range 
azimuth 
elevat i on 
round  number 
radius  of  the  earth 
range 

type  of  rocket 
station  altitude 
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SHEAR 

SIG  (4,200) 
el  t  1 
2 

3 

4 

SITE  (3) 
SLAT 
SLAT  2 
SLAT  4 
SLONG 

SMOT  (10,5) 
elt  1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

SM4 

SM5 

SS  (5,6) 
elt  1 
2 

3 

4 

5 

STAT 
SZ  (11) 

TEMP  (200) 

TFF 

TIE 

TIEP 

TIER 

TIET 

TIME  (200) 

TIMEL 

TIMEZ 

TIMM 

TM 

TP  (4,200) 
elt  1 
2 

3 

4 

TPCOR  (200) 
TT 


wind  shear 

array  containing  data  for  significant  lev 

corrected  temperature  °K 

uncorrected  temperature  °C 

pressure 

alii  tilde 

launch  site 

station  latitude 

SLAT  multiplied  by  2 

SLAT  multiplied  by  4 

station  longitude 

array  containing  smoothed  radar  data 

time 

range 

azimuth 

elevat ion 

X 

Y 
Z 

X  velocity 

Y  velocity 
Z  velocity 

1-eharacter  field  containing  sign  >>t 
N-S  wind  c  imponent 

1-character  field  containing  sign  of 
E-W  wind  component 

array  containing  input  base- level  points 

al ti tude 

pressure 

raob  temperature 

rocket  temperature' 

temporary  storage- 

station  number 

array  containing  altitudes  of  mnstant- 

pressure  levels 

array  containing  temperatures 

first  time  read  from  radar  input  tape 

base- level  rocket  temperature  °G 

base- level  pressure 

base- level  raob  temperature 

base-level  rocket  temperature  °K 

array  containing  time's 

time  of  launch  (GMT) 

lift-off  time  (sec) 

time  in  minutes  and  seconds 

time  of  raw  data  point 

array  containing  input  temperature  data 

time 

temperature 
a  1 1 i tude 

temperature  correct i on 

array  containing  temperature  corrections 
time  of  wind  data  point 
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TT1 

VD 

VD1 

VS 

VS1 

VT 

VT1 

VT2 

VT3 

VX 

VX1 

VY 

VY1 

VZ 

VZP  (200) 
WGT  (59) 

WSSITE  (3) 
WSSTAT 
WW  (10) 

XAC 

YAC 

YK 

ZAC 

ZF 

ZFIRST 

ZLAST 

ZZ 

ZZ1 

MKN 

B 

F9 

F10 

SM10 

SM7 

SM8 

SM9 


time  of  previous  point 
wind  direction 

wind  direction  of  previous  point 
windspeed  . 
windspeed  of'  previous  point 
total  velocity 

total  velocity  of  previous  point 

VT1  multiplied  by  1.3 

VT1  multiplied  by  .7 

corrected  E-W  wind  component 

E-W  wind  component  of  previous  point 

corrected  N-S  wind  component 

N-S  wind  component  of  previous  point 

vertical  velocity 

array  containing  fall  velocity 

weighting  factors  used  in  the  117-point 

raw  data  filter 

alpha  field  containing  SMR,  WSMR,  NM 

station  number  for  WSMR 

weighting  factors  used  in  the  10-point 

filter  used  to  smooth  accelerations 

smoothed  X  acceleration 

smoothed  Y  acceleration 

year 

smoothed  Z  acceleration 
altitude  in  feet 

altitude  of  first  valid  wind  data  point 
altitude  of  last  valid  wind  data  point 
altitude  (meters)  corrected  for  earth's 
curvature 

altitude  of  previous  point 


temporary  storage 

pressure  divided  by  power  of  10 

density  divided  by  power  of  10 

sign  of  exponent  of  10  used  with  density 

sign  of  temperature 

sign  of  temperature  correction 

sign  of  exponent  of  10  used  with  pressure 


WINDAVE 

CON 

D1V 

JL 

KK 

KJ 

NOM 


factor  to  convert  altitude  to  meters 
averaging  interval  ize 
=  0  if  wind  input  is  radar  tape 
equal  2  if  altitude,  to  be  used  are  in 
meters,  3  for  feet 

equal  0  until  first  point  has  been  processed 
number  of  averaged  points 
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NZ 

number  of  consecutive  bad  wind  data  point 

TPREV 

time  of  previous  point 

TT 

time 

VD 

wind  direction 

VS 

windspeed 

VX 

E-W  wind  component 

VY 

N-S  wind  component 

VZ 

fall  velocity 

ZPREV 

altitude  of  previous  point 

ZZ  (2) 

array  containing  data  point  altitude 

elt  1 

meters 

2 

feet 

RTDATA 

BUFFA  (3241) 

input  data  buffer 

IA  (32) 

array  containing  azimuth  data  by  subchannel 

ID  (32) 

array  containing  rate  of  change  of  range 
(R  dot)  data  by  subchannel 

IE  (32) 

array  containing  elevation  data  by  suluhanne 

I  EOF 

end  of  file  indicator 

I  FIRST 

equals  0  until  initialization  is  complete 

IMTDAT 

array  used  when  reading  met  data 

INBUF  (3241) 

input  data  buffer 

I  READ 

end  of  data  block  indicator 

1TQ  (32) 

array  containing  tracking  mode  data  by 
subchanne l 

1TYPE 

data  type  indicator 

KK  (6) 

status  array  used  by  tape  read  routine 

KNT 

number  of  points  of  data  processed  in  data 

block 

T 

time  of  data  point 

T1 

time  of  met  data  point 

TPREAD 

AZM  (32) 

array  containing  azimuth  data  by  subchannel 

A1 

azimuth  value  used  in  2-point  average 

A2 

azimuth  value  used  in  2-point  average 

ELV  (32) 

array  containing  elevation  data  by  subchanne 

ENDFI1 

end  of  file  indi cator 

F.1 

elevation  value  used  in  2-point  average- 

E2 

elevation  value  used  in  2-point  average- 

ID  (32) 

array  containing  radar  identification  by 
subchannel 

I  PAR 

parity  error  indicator 

IRDT  (32) 

array  containing  rate  of  change  of  range 
(R  dot)  by  subchannel 

[TQ  (32) 

array  containing  radar  tracking  mode  by 
subchannel 

J1 

index  for  data  for  first  radar  processed 

J2 

index  for  data  for  second  radar  processed 
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RAN  (32) 
R 1 
R2 

TALT 
T 1 

THERMO 

JJ 

KKK 


array  containing  range  data  by  subchannel 
range  value  used  in  2-point  average 
range  value  used  in  2-point  average 
time  ot  data  point 

time  value  used  in  2-point  average 


index  of  element  in  array  ARC  which 
corresponds  to  jth  element  in  ALT  array 
index  of  lirst  element  for  which  pressures 
are  computed 


D 

DZ  (46,3) 

RATIO 
2  (46) 


OZONE 

01 

02 

ZX 

ZY 


interpolated  value  of  DZ 

array  containing  coefficients  used  in 

Thermo  correction 

constant  used  in  interpolation 

array  containing  altitudes  corresponding 

to  coefficients  AZ ,  BZ ,  CZ,  and  DZ 


d(ln  l >)/dz 
d( In  T)/dz 

change  in  altitude  from  following  point 
change  in  altitude  from  previous  point 


VIII.  TEMPERATURE  REDUCTION  COMPARISON 

Since  this  report  describes  a  system  that  will  perform  automated 
processing  of  the  temperature  data  from  a  rocketsonde,  an  indication 
of  the  accuracy  of  such  processing  will  now  be  made. 

Because  the  temperature  data  from  the  automated  system  merely  re¬ 
places  the  manual  data  with  no  further  changes  in  the  computer  process¬ 
ing  of  the  rocketsonde,  only  comparisons  of  temperatures  from  both 
systems  need  be  made. 

As  an  indication  of  accuracy  between  the  two  systems,  three  sepa¬ 
rate  rocket  flights  were  processed  by  the  automated  system  and  the 
temperature  data  compared  with  manual  data  in  graphical  plots.  These 
results  are  shown  in  Figures  10,  11  and  12. 
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In  comparing  the  two  plots  in  each  figure,  there  is  good  agree¬ 


ment  in  some  areas  and  not  so  good  agreement  in  other  areas.  Since 
the  criterion  specified  by  the  IR1G  data  reduction  requires  that 
the  straight  line  segments  he  only  within  2°  of  the  actual  curve — there 
curve--there  is  some  room  for  difference.  For  example  in  Figure 
Is  some  room  for  difference.  For  example  in  figure  10,  the  point  at 
2150  s  and  -58°  yields  a  difference  of  3.5.  If  the  actual  (strip  chart) 
temperature  were  -59.75°,  then  the  differences  (manual  and  computer 
verses  actual)  would  both  be  1.75°,  well  with  the  2nd  degree  tolerance. 

In  a  further  analysis  for  comparison,  the  actual  raw  data 
plots  (Q5  recording)  were  compared  with  both  the  manual  and  auto¬ 
mated  output  of  points  where  the  differences  were  large.  In  all 
cases  except  one,  both  data  were  within  the  2°  limit.  The  one 
exception  was  the  manual  reduction,  which  was  in  error  due  to  the 
incorrect  reading  of  the  raw  data  chart. 

These  comparisons  indicate  that  the  automated  reduction  method 
is  as  good  or  better  than  the  manual  technique. 

IX.  DESCRIPTION  OF  DATA  CARDS  FOR  PROGRAM 
A.  IDENTIFICATION  CARD 

The  identification  card  is  included  in  all  runs  and  contains 
general  information  about  the  type  of  round  and  the  launch  location. 

Col  No  Format 
01-18  3A6 

19-22  A4 

23-24  A2 

25-26  21 

27-28  A2 

29-32  A4 


Rocket  Type  (Super  Loki ,  etc.) 

Round  number  for  the  type  of  rocket 
fired  during  the  year  from  that  site 
Last  two  digits  of  year 
Two-digit  number  representing  month 
Two-digit  number  for  day  of  month 
GMT  time  of  launch  rounded  to  nearest 
minute  (leave  remaining  fields  blank  if 
launch  site  is  SMR,  WSMR 

Name  of  launch  site,  e.g.  ,  Green  River,  Utah 
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33-50  3A6 


31-3S 

Ar> 

MRN  statu  n  number  (e.g 
River) 

atj-t)4 

F9 . 4 

Station  latitude  (e.g., 
River) 

co-73 

F9 . 4 

Station  longitude  ( e . g . 
R i ve  r  ) 

74-78 

1-7.2 

Station  altitude  meters 
Green  River) 

,  724*7  tur  lir'-t-n 
38.933  fur  Green 
1 1 0 . Oh  j  inr  '  i  r t •  t ■  1 1 
MSI.  (e.g.,  1  50,s  I  •  ■ 


B.  OPTIONS  CARD 


Ttif  options  card  is  included  in  nil  runs  .un!  is  used  lor  n'lri- 
ting  the  options  for  each  run. 


Col  No  Format 

1  1  1 

I  OP  (  1 ) 


:i  1 1 


IUP(2) 


5  1 1 

lOPCi) 


7  1  1 


I0P(4) 

9  1 1 

!0P(5) 


Wind  input  type 

=  (0  or  blank  i  -  no  wind  input,  temper. i- 
tures  only 

-  1  UN l LOG  radar  tape  used  to  compute  wands 
=  2  DR  format  radar  tape  used  to  compute  wm 
=  3  winds  input  on  cards 

Temperature  input  type 

=  (0  or  blank)  -  no  temperature  input  -  wind 
only  (no  temp  sense) 

-  I  time  vs.  temp  input  on  carils  (used  only 
with  winds) 

=  2  altitude  (km)  vs.  temp  input  on  cants  ( it 
wind  input) 

=  3  altitude  (ft)  vs.  temp  input  on  cards  <  :i 

wind  input) 

=  4  no  temperature  input  (sensor  failure! 

=  5  temperature  vs.  time  from  digitized  tele 
me  t  ry 

Temperature  correction  type 
=  (0  or  blank)  -  no  temp  correction  used 
=  1  Arcasonde  temp  correction  used 
=  2  Loki  Datasonde  temp  correction  used  ((or 
runs  with  wind  and  temperature,  for  tempera¬ 
ture-only  runs,  no  temperature  correction  is 
computed  but  one  may  be  included  on  the  tem¬ 
perature  data  cards). 

=  3  Krumins'  correction 
English  units  option 

=  (0  or  blank)  -  no  additional  output 
f  (0  or  blank)  -  additional  listing  and  out¬ 
put  cards  containing  data  in  English  units 
for  1000  ft  levels 
Ozonesonde  parameters  option 
=  (0  or  blank)  -  no  additional  output 
f  (0  or  blank)  -  data  computed  for  . 9  km 
and  1  kfl  instead  of  standard  1  km  and 
5  kft,  ozonesonde  parameters  computed  bu 
all  levels 
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!  1  II 

I  UP ( b ) 
!  5  11 


1 0P( 7 ) 


11 


I  1 


I  OP ( 8 ) 


17  II 

I OP (9) 


19  I  1 

10P(10) 


21  I  1 

IOP(ll) 


Ignore  wind  data  i oust ra i nt s  opt  i  .n 
=  (0  or  blank)  -  stamlar'l 

t  (0  or  h  1 . : :  i  k  )  -  imisi  r.iints  on  fall  U'l"(  it\ 
art<l  point-to-j)oint  wiini  variation  ignored 
Corrected  winds  used  on  card  input 
=  (0  or  blank)  -  assumes  input  wind  data  is 
uncorrected 

t  (0  or  blank)  -  assumes  input  wind  data  is 
corrected 

(This  option  is  used  only  when  winds  are 
input  on  cards.  The  only  difference  in  out¬ 
put  is  the  columns  in  which  the  wind  is  put 
in  the  MRN  formatted  output). 

Ignore  base-level  criterion  option 
=  (0  or  blank)  -  if  no  input  base-level 
point  has  agreement  of  2.5°  between 
raob  and  rocket  temperatures,  no  thermodyna¬ 
mic  data  is  computed.  11  more  than  one  point 
meets  the  criterion,  the  one  closest  to  25  km 
is  selected. 

t  (0  or  blank)  -  if  one  or  more  points  meet 
the  above  criterion,  the  one  closest  to  25 
km  is  selected. 

Continue  temperature  data  option 
t  (0  or  blank)  -  if  a  layer  of  missing  tem¬ 
perature  data  is  greater  than  5  km  -  tem¬ 
peratures  are  considered  missing  for  remain¬ 
der  of  f 1 i ght . 

f  (0  or  blank)  -  temperature  computations 
are  continued  regardless  of  existence  of 
missing  data 

Continue  thermodynamic  data  option 
=  (0  or  blank)  -  if  a  layer  of  missing 
temperature  data  is  greater  than  ]  km, 
pressure,  density,  etc.  .ire  not  computed 
for  remainder  of  flight 

f  (0  or  blank)  -  thermodynamic  computations 
are  continued  regardless  of  existence  o| 
missing  data. 

Suppress  printing  1/sec  wind  data 
=  (0  or  blank)  -  standard 

f  (0  or  blank^  -  1/sec  wind  data  points  are 
not  printed.  The  summarized  output  remains 
unchanged . 


C.  DIGITAL  RADAR  TAPE  PARAMETERS 

This  card  is  included  only  when  winds  are  to  be  computed  from  radar 
data  input  on  a  digital  tape  (i.e.,  UNILOG  tape  or  DR  formatted  tape). 
Col  No.  Format 
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1-3 

13 

Radar  numl  '-r 

5-6 

12 

Radar  subchannel 

(  1  eave  b  1  arik  for 

DR  tap* 

8-10 

) 

Radar  number  for 

2nd  Rada  r  (  l  t-avt 

blank  1 

DR  tape  or  for  pr 

ocessing  only  one 

-  r.ida  r  ' 

12-13 

12 

Radar  subchannel 

fur  2nd  radar 

15-23 

K9 . 3 

Lift-off  time  in 
tape) 

seconds  ( ! caVf  b 1 

ank  Yu 

25 

1  1 

Number  of  reels  of  input  tape 

D.  TEMPERATURE  CALIBRATION  DATA 

These  cards  are  included  only  when  temperature  input  option  and 
temperature  vs.  time  t rom  digitized  telemetry,  is  selected.  The 
data  for  these  cards  is  taken  from  the  calibration  chart  data. 

1 .  Ca  rds  1 ,  2 ,  and  3 

Col  1  thru  70  7K10.4  format  -  the  21  values  located  in  the 
2nd  column  of  the  calibration  sheet  that  shows  the  instru¬ 
ment  number  and  thermistor  number  for  the  datasunde.  See 
Table  4. 

2.  Card  4 

Col  1  thru  50  5F10.3  format  -  the  5  values  located  in  the  last 
row  of  the  sheet  used  for  cards  1,  2  and  3  above.  See  lable  a. 

3.  Card  5 

Cot  1  thru  30  3F10.3  format  -  the  3  values  of  temperature 
correspond i ng  to  the  ratios  of  0.3,  0.5,  and  0.7,  respectively, 
from  the  ratio  vs.  temperature  data. 

4.  Card  6 

Col  1  thru  10  F10.1  format  -  this  value  is  taken  I rom  the 
Q-5  plot  and  represents  the  ordinate  corresponding  to  r«-‘Yreme 
during  the  T-2  minute  prelaunch  plot.  This  value  is  not  neces¬ 
sary  for  temperature  calculations  but  0  >es  provide  a  cross 
reference  value  to  the  Q-5  recording. 

F. .  WIND  DATA  CARDS 

These  cards  are  used  to  input  manually  reduced  winds  and  art- 
omitted  if  a  radar  tape  is  used  or  for  a  Lemperatures-only  run. 

There  is  no  program  limitation  on  the  number  of  wind  data  cards 
which  may  he  included  in  any  run. 


Col  No 

Format 

1 

1 1 

=  1  card  id  for  wind  cards 

2-10 

F9 . 6 

Time  from  lift-off  (min.  sec.) 

11-22 

F10.2 

A 1 1  i  t  ude  km  MSI. 

21-30 

F10.0 

Wind  direction  -  degrees  1 rom  north 

31-40 

F10.0 

Windspeed  -  meters/sec 
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F.  BASK- l.KVF.L  POINT  CARDS 

Those  cards  art'  used  to  input  a  base- 1  eve  1  pressure  arid  corres¬ 
ponding  altitude  from  which  thermodynamic  data  (pressure,  density) 
is  computed  using  hydrostatic  equations  and  rocket  temperatures. 

Zero  to  5  base  level  points  are  input  with  each  run  with  temperatures 
and  the  best  one  is  selected  as  described  for  the  8th  option  on  the 
option  card.  If  no  base-level  point  is  input  or  it  all  art-  rejected 
because  of  temperature  disagreement,  then  pressures  and  densities 
arr>  not  computed. 

Col  No  Format 

111  -2  card  id  for  lise-level  point 

2-10  F9.2  =  Geometric  altitude  -  km  MSI.  lit  MSI. 

are  used  if  temperature  altitudes  are 
i nput  i n  ft  ) 

11-20  F10.3  Raob  pressure  -  millibars 

21-30  F10.3  Raob  temperature  -  °C 


G.  TEMPERATURE  DATA  CARDS 

These  cards  are  used  to  input  temperature  on  all  runs  that 
include  temperatures  when  the  temperature  vs  time  -  telemetry  option 
has  not  been  selected.  As  many  as  200  temperature  data  cards  may  bl¬ 
ind  uded  in  any  run. 

Col  No  Format 


1 

11 

=  U  card  id  for  temperature  cards 

2-10 

F9.2 

Time  from  lift-off  (min.  sec.) 

11-20 

F10.3 

Temperature  -  °C 

(Rest  of 

card  is  blank  except 

for  tempera t ures-on 1 y  runs) 

21-30 

F10.3 

Geometric  altitude  -  km  or  It  MSI. 

-  (see  2nd  option  on  options  card) 

31-40 

FI  0.3 

Temperature  correction  -  °C 
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x.  examples 

A.  I'EMI’EKATUKK  VS.  T I ME  ERuM  UK.ITIXEU  1 I.  I.EMETRY 

An  examp  1 1‘  of  how  the  pro)' ram  can  In-  rim  from  a  lrnmii.il  ■■■  i  i  ! 
now  l>c  given. 

The  Table  4  in*  1  ow  shows  most  ol  ttn-  c  a  I  i  brut  i  on  data  ttiat  n 


needed  ‘nr  computing  temperature  versus  time  from  telemetry 

( I0I’(2 )  =  5)  . 

The  data  within  ttie  rectangle  (Column  2)  of  Table  4  is  input  ■ 
41  thru  61,  and  the*  data  within  the  rectangle  (last  row)  of  1  ,j  1  >  1  < -  a 


is  input  on  71.  See  next  page. 

The  remainder  of  the  data  is  as  described  in  Section  Y 

INSTRUMENT  NO  21540  THERMISTOR  NO  77-10-10212 
2  I 


1.0 

0.9915 

5572.0 

5  b  20 . 0 

5.0 

0.9577 

i  5572.0 

5818.. 

10.0 

0.9190 

5572.0 

6063.0 

20.0 

0.8509 

5572.0 

6548.0 

30 . 0 

0.7935 

5572.0 

7022.0 

40.0 

0.7440 

5572.0 

7489.0 

50.0 

0.7006 

5572.0 

7953.0 

60.0 

0.6623 

5572.0 

8413.0 

70.0 

0.6299 

5572.0 

8846.0 

80.0 

0.6000 

5572.0 

9286 . 0 

90.0 

0.57.36 

5571.0 

9713.0 

100.0 

0.5506 

5571.0 

10118.0 

150.0 

l  0.4594 

5571.0 

12126.0 

200 . 0 

0.3986 

5571.0 

13976.0 

300.0 

0.3222 

5571.0 

17291.0 

400.0 

0.2766 

5571.0 

20139.0 

600.0 

0.2234 

5571.0 

24938.0 

800 . 0 

0.1941 

5571.0 

28708.0 

1000.0 

0.  1757 

5572.0 

31722.0 

1500.0  i 

0.  1498 

5572.0 

.37191.0 

2000.0  | 

,  0.1359 

5572.0 

40986.0 

.-65 -0 

-50.0 

-25.0 

0.0 

15.0 


1333.000  479.600 


109.400 


31.000  7.107 


TABLE  4  THERMISTOR  ANT)  INSTRUMENT  CALIBRATION  DATA 
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TERMINAL  COMMANDS 


NOTE  ALL  LINES  THAT  HAVE  >  IN  THEM 
ARE  TERMINATED  WITH  A  RETURN  ON  TERMINAL 

>@USE  A. ,8300-DATAFL. 

READY 
>PF.D ,  I  A. 

CASE  UPPER  ASSUMED 
ED  16R1W1  -  (DAY  -  TIME) 

INPUT 
II  :> 

EDIT 

0:>TAB  /  This  establishes  the  tab  character 
0 : >SET  19 
0 ;  > 

INPUT 

1 1 : > LOK I  PARACHUTE/00137901221900 
21  :> 

EDIT 

1 : >SET  3,3,7,17,19,21 
1  :> 

INPUT 

21 : >  1/5/2/  1/ 1 / 1 
31  :> 

EDIT 
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2: >SET  5, 8, 15,25 
2 :  > 

INPUT 

31  :  >112/1/296/68400.0/1 
41  :> 

EDIT 

3:  'SET  11,21,  11,41,51,61 
.1 :  > 

INPUT 

41  :  >.99 15/. 9577/. 9 190/. 8509/. 7935/. 7440/. 700t. 
5 1 : > . 6623/ . 6299/ . 6000/ . 5736/ . 5506/ .4594/ . 3986 
61 :>.3222/ .2766/. 2234/. 1941/. 1757/ . 1498/  .  1359 
71 : >1333-/479.6/ 109.4/31 ./7. 107 
81  :>-44.5/-27.6/-10. 

91 : >90.2 

101: >2  20.939/50. / - 5 5 . 6 
111 :>2  24. 021/30. /-49. 2 
1 2 1 : >2  26.953/20. /-46. 9 
131 : >2  31 .577/10. / -37 . 8 

141  :> 

EDIT 
1 3 : >LNP ! 
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l.-LOKI  PARACHUTE  00587901221900 


2:1521  11 

1:111  1  196  68400.0  1 


4:  .9915 

.9577 

.9190 

.8509 

.7935 

.  7440 

.  7000 

5:  .  6623 

.6299 

.6000 

.5736 

.5506 

.4594 

.  3986 

6:  .3222 

.2766 

.2234 

.  1941 

.  1757 

.  1498 

.  1359 

7: 1333. 

479.6 

109.4 

31.0 

7.107 

8: -44.5 

-27.6 

-10. 

9:90.2 

10:2  20.939 

50. 

-55.6 

11:2  24.021 

30. 

-49.2 

12:2  26.953 

20. 

-46.9 

13:2  31.577 

10. 

-37.8 

TABLE  5 

Card 

Images  for 

ROCKET 

Program 

0 :  >F.XIT 

LINES:  13  FI  ELD AT A 

(ANY  ERRORS  MADE  ABOVE  CAN  BE  CORRECTED  BY  USING  THE  OTHER  EDIT  COMMANDS') 
>(<iUSE  B.  ,8300'"R0CKET. 

READY 

>PASG,TZ  2.,U,ACB788  (Assign  unilog  tape) 

>@ASG,ZA  PUNCHFL. 

>@ASG,ZA  PRINTFL. 


>('JBRKF 


TNCHS/PUNCHFL 


>(uBRKPT  PR1NT$/PRI.NTFL 

>(>'HDG  (NAME  NUMBER  FOR  OUTPUT) 

>(->XQT  B .  NETROC 

>foAI)D  A. 

(There  will  be  a  wait  period  here  until  a  >  is  printed) 

>PBRKPT  PUNCH  S 
>PBRKPT  PR 1NT$ 

The  results  of  the  output  PR1NTFL  ran  he  sent  to  any  of  tin- 
high  speed  printers  by  use  of  the  SYM  control  statement  or  the  user 
ran  examine  the  output  with  the  El)  processor.  For  this  example-  the 
results  of  the  computations  up  to  where  the  wind  data  is  printed  is 
shown  on  the  next  pages. 
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OPTIONS  SELECTED 

WIND  INPUT -UN I  LOG  TAPE 

TEMPERATURE  VS  TIME-TELEMETRY 

HENRY  ( IR1G )  CORRECTION  FOR  LOKI  DATASONDE  TEMPERATURES  USED 
CONTINUE  THERMODYNAMIC  COMPUTATIONS  EVEN  IF  LAYER  OF  TEMP  DATA  MISSING 
IS  GREATER  THAN  3 KM 
SUPPRESS  PRINTING  1/SEC  WIND  DATA 

LOKI  PARACHUTE  NUMBER  0013  LAUNCHED  22  JAN.  79,  1900  from  SMR ,  WSMR,  NM 
STATION  LATITUDE  32.467  LONGITUDE  106.417  ALTITUDE  1215.910  METERS 
RADAR  113 

GEOMETRIC  ALTITUDES 


***LOKI  CALIBRATION  TAPE  DATA*** 
*Rat i os* 

.9915 
.9577 
.9190 
.8509 
.7935 
.7440 
.7006 
.6623 
.6299 
.6000 
.5736 
.5506 
.4594 
.3986 
.3222 
.2766 
.2234 
.  1941 
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.  1  7  r>  7 
.  1498 
.  1.189 

•■CALIBRATION’  VALUES* 

1111. 0000  479.6000  109.4000  11.0000  7.1070 

•'•CALIBRATION  CHECK* 


RATIO 

COMP 

CHART 

D I  FF 

.3 

-44.7 

-44.5 

-  .  1 

.  5 

-27.4 

-27.6 

_  > 

.7 

-10.0 

-10.0 

.0 

DATA  FOUND  IN  SUBCHANNELS 
1 
9 
19 

FIRST  TIME  ENCOUNTERED  ON  RADAR  TAPE  .068  SEC 


LOKI  PARACHUTE  NUMBER  0013  LAUNCHED  22  JAN  79,  1900  FROM  SMR,  WSMR,  NM 
STATION  LATITUDE  32.467  LONGITUDE  106. 4 1 7  ALTITUDE  1215.910  METERS 
RADAR  113 

GEOMETRIC  ALTITUDES 

•'•'SIGNIFICANT  LEVEL  TEMPERATURES* 


TIME 


SEC 

MIN:  SEC 

TEMP 

RATIO 

ORDINATE 

123 

2:03 

31.71 

.9299 

83.7 

127 

2:07 

33.99 

.9355 

84.2 

156 

2:36 

.32 

.7905 

71.0 

166 

2:46 

-3.21 

.7786 

70.0 

175 

2:55 

-8.79 

.7122 

64.0 

185 

3:05 

-8.16 

.7184 

64.6 

195 

3:  15 

-14.30 

.6554 

58.9 

213 

3:33 

-15.76 

.6392 

57.4 

2.32 

3:52 

-12.32 

.6761 

60.7 

259 

4:  19 

-17.37 

.6210 

55.7 
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2  79 

4 :  39 

-11.08 

.6891 

61  .9 

293 

4:53 

-9.30 

.7072 

6  3 . 6 

315 

5:15 

-15.94 

.6373 

57.3 

330 

5  :  10 

-13.20 

.  6669 

60.0 

367 

6:07 

-21  .  10 

.  5769 

51.9 

401 

6:41 

-20.86 

.5802 

52.2 

4 1 4 

6 : 54 

-15.86 

.6387 

57.5 

443 

7:23 

-12.04 

.6795 

61.2 

465 

7:45 

-12. 14 

.6787 

61 . 1 

475 

7:55 

-17.46 

.6202 

55 . 8 

484 

8:04 

-31.76 

.4466 

40.2 

489 

8:09 

-34.62 

.4111 

37.0 

516 

8 :  36 

-36.23 

.  3920 

35 . 3 

535 

8:55 

-31.42 

.4500 

40.5 

567 

9:27 

-37.68 

.3767 

33.9 

596 

9:50 

-35.36 

.4024 

36.3 

623 

10:23 

-39.13 

.  3606 

32 . 5 

632 

10:32 

-37.77 

.3753 

33.9 

674 

11:14 

-43.11 

.3166 

28.6 

707 

11:47 

-41.43 

.  3345 

30.2 

744 

12:24 

-45.44 

.2930 

26.4 

804 

13:24 

-42.51 

.3224 

29.1 

836 

13:56 

-46.05 

.2873 

25.9 

911 

15:11 

-48.35 

.2662 

24.1 

974 

16:14 

-43.31 

.3142 

28.4 

1058 

17:38 

-39.67 

.3542 

32.0 

1 1 35 

18:55 

-44.24 

.3056 

27.7 

1176 

19:36 

-41.62 

.  3325 

30.  1 

1227 

20:27 

-43.62 

.3117 

28.2 

1278 

21:18 

-41.27 

.3369 

30 . 5 

1333 

22:13 

-46.45 

.2836 

25.7 

1485 

24:45 

-45.95 

.2882 

26.1 

1600 

26:40 

-51.81 

.2376 

21.5 

1783 

29:43 

-49.20 

.2595 

23.5 

2223 

37:03 

-55.99 

.2059 

18.6 

2403 

40.03 

-53.46 

.2234 

20.3 
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B.  TEMPERATURE  VS.  TIMF.  FROM  CARDS 

On  page  60  is  .1  listing  of  a  file  called  8300'-vL0K  I  .  CONTROL . 

This  file  contains  information  for: 

1)  Control  of  operations  -  Lines  1-8  anil  Lines  a4-45 
1)  Data  for  IDENTIFICATION  CARD  (page  46)  -  Line  4 

3)  Data  h  r  OPTIONS  CARD  (page  47)  -  Line  10 

4)  Data  for  DIGITAL  RADAR  TAPE  PARAMETERS  (page  48)  Line  - 

8)  Data  for  BASE-LEVEL  Point  Cards  (page  50)  -  Lines  12-15 

6)  Data  for  TEMPERATURE  DATA  CARDS  (page  50)  -  Lines  lo-43 


Prior  to  execution  of  the  control  file,  a  phone  call  to  the  1 
library  (678-3173)  requesting  that  tape  E206  he  taken  to  system  1! 
should  he  made.  The  execution  is  made  by  typing  (the  '  '  is  the  ; 
premier ) 

•KADI)  8JOO*LOK I.  CONTROL 

Whiui  the  system  responds  with  the  file  PRINT,  can  he  examinei 

with  the  El)  processor  or  -SVMmed  to  a  printer. 

Re  low  is  an  explanation  of  the  listing  shown  on  the  pages 
nidi ca ted : 


age 


LI. STING 


til 
62 
h  1 
64 
6.5 
66 

67 

68 
6<) 

70 

71 


Options  selected 

tV,,  per  second  wind  (would  continue  until  end  of  I 
Thermodynamic  data  for  input  temperatures 
Thermodynamic  data  for  kilometer  Levels 
R0C0B  message 

Plot  of  Wind  (x,y)  and  Temp  vs.  Altitude 
MRN  30  cards 

MRN  30  cards  (continued) 

Thermodynamic  data  for  5000  feet  levels 
Thermodynamic  data  for  mandatory  levels 
MRN  40  cards 


1 1 


ape 


iv  St  t  'III 


i  ght  ) 
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«!  * 


10 
1  l 
1  2 
i  3 
l  * 
lb 
1  6 
17 

1  b 
J  V 
20 

2  I 
22 
23 
2*4 

25 

26 
'  2> 

2b 

29 

30 

3  1 
32 


b  300 •  L  0 *  l  CONlfiL 


19  tU  6  1 


PfrlNT. 

bbKM'l  PH1».  I  *  /  H  N  1  NT 


„  A  b  U  ,  C  P  KuM.ru 
Wtt^APT  KU'.C  «S  /  PONC  h 
•LNAbCt  ® 3Uu  • KOC  K  £  1  • 
*AbO,A  tt3GO*H0CA£T# 

•  Abu,  If  2  ,  u  ,  i 2  O  6 

•  xwl  6  JUO«  W  uc  A  t  T  •  HE  T  a  OC. 
’  L  OH"  FTk  aChu  T  t 

1  i  3 

1  M  uz 

22  l  .61 1 
22^-060 
22*»  .  6bV 


oL  v2o 1 UV25 I  6  30 


22b  .  2 W 
1  2.39 

1  3  .  0  1 

1  3.3b 
1  3.16 

1  1.2b 


1  1.16 
1  b  .  06 
1  b  .  9  tt 
1  6.0V 
1  7,39 

9  7.52 

1 
1 


tt  .  33 
v  .  b  V 

lo.  37 
11.30 
I3.|2 
1  b.  2  9 


59100.151  1 

9  9.6 
36.  V 
26.5 


25.3 
-26.2 
-15.  ? 
-19.7 

-  I  7.  J 
-08.  I 
-09  .  b 
-07.1 

—  10.6 

-  U  9 . 9 

-  I  3.  b 

_-l0.2_ 

-09.3 
-2  1  .b 

-  26.6 
-32.6 
-32.3 
-36.9 


3T 

- T 

T5TJ7 

39 

9 

i  a .ob 

JS 

9 

21.36 

-99 . 3 

36 

9 

22. bV 

-98 . 9 

37 

9 

27. bb 

-b2.i 

38 

9 

30,32 

-S2.» 

3  9 

9 

J  2  •  J  f 

"-55 

9  U 

9 

3».U/ 

-bb.  a 

91 

9 

36  .U2 

-be .  ii 

92 

9 

36, b3 

-57.9 

93 

9 

36, blf 

-56.8 

99 

•  BKAPT  PNlttls 

b  7  •  e 
5bl9 
b2. 7 


■  b  2  •  2 
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XI.  TAPK  UNIT  ASSIGNMENTS 
INPUT  TAPES 

1.  UN  I  LOG  Tape  -  radar  data  recorded  on  a  digital  tape  using  a  real¬ 
time  log  routine  on  the  UN  I VAC  1108.  The  UNI  LOG  Tape  is  assigned  to 
Fortran  logical  unit  2  and  is  read  hy  the  system  subroutine  FLDT'AP. 

Each  physical  record  contains  data  for  a  one-second  interval  and 
begins  with  the  first  message  after  the  whole  second.  The  data  rate 
is  20  messages  per  second.  The  record  length  is  variable. 

Each  physical  record  contains  a  control  word  and  20  logical 
records.  The  control  word  should  be  ignored.  Each  logical  record 
contains  a  word  count  and  a  time  word,  which  are  followed  by  five 
words  of  information  for  each  active  input.  The  maximum  number  <>t 
inputs  is  32.  Therefore,  the  maximum  size  of  a  physical  record  is 
3241  words  (1  +  (.20"  (2  +  (32v'-5))  )).  If  no  inputs  are  active,  each 
logical  record  will  contain  only  a  word  count  (equal  to  one)  and  a 
time  word.  The  minimum  size  of  a  physical  record  is,  there* fore,  41 
words  (1  +  (20" (2)  )). 

The  time  word  is  floating-point  range  time  in  milliseconds.  The 
five  data  words  for  each  input  contain  information  from  the  standard 
serial  data  format.  The  IOMUX  input  subchannel  number  appears  in  the 
second  sixth-word  (second  most  significant  character)  in  the  second  data 
word  of  each  five-word  input  data  block.  This  number  has  a  value  from 
zero  to  31.  An  example  of  a  physical  record  containing  radar  data  is 
on  the  following  page. 
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UN  I  LOG  TAPK  KuKMAT 


CONTROL  WORD 
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1 

35  33  32 
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29  24 
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pRun  %  (  TGT 
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Tsyvc  i  nr 

1000  It  m  1  , 
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o  ! 

DATA 
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1  1  s/c 
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LSB  1  1 
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o  ; 
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Tmsb 
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WORD  4 
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ELEVATION  1 

LOGICAL 

23 
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RECORD 

DATA 

WORD  5 

L_ 
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RANGE  RATE  LSB  I 

35  33  32 

30 

29  24 

23  16  15  9 

8  5  0 

PHYSICAL 

DATA 

WORD  1 

[  RUN  if  ]~  TGT 

»  T 

C  C 

!  SYNC'  1  ID 

ToooTT  Ml 

RECORD 

29 

24 

23 
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DATA 

WORD  2 

1.  I  S/C  if* 
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RANGE 

LSB] 

23 
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DATA 

WORD  3 

[ _ 

Tmsb 

AZIMUTH 
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23 
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DATA 

WORD  4 

r~ 

Imsb 

ELEVAT 1  ON 
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23 
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DATA 

WORD  3 
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1 

RANGE  RATE 

LSB 

19  more 

logical  records 

TGT  il  =  Target  Number 


CC  =  Confidence  Counter 
II)  =  Site  Identification 
TM  =  Tracking  Mode 

*  Subchannel  Numbers  will  be  in  ascending  order  (0  -  31). 
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I  lie  rang*'  is  1  ii  yards  and  the  azimuth  and  e  I  ev.i  l  i  ori  are  i  n  units 
equivalent  to  4.79  57x10  ^  radians. 

2.  1  j R  format  tape  -  radar  data  recorded  on  a  digital  tape  gen- 


eraled  by  pr 

oi ess i ng  the  rada r  fie 

Id  tape  DR  formatted  tapes 

ass i gned 

to 

Fort  ran  1 og i c  a  1  un i t  1 

These  tapes  are  BCD  and  i 

into  nti.i  t 

i  on 

used  hy  MF.TH0C  in  the 

following  t  o  rmu t . 

‘!h. t  carte 

r  No 

Format 

11-19 

F9 .  ii 

Time  (seconds  from  lift-off 

a  1  -  V) 

F9 .  1 

Slant  range  from  radar  (yar 

1,1-99 

F9 . 2 

Azimuth  from  radar  (mils) 

71-79 

F9 . 2 

Elev.it  ion  from  radar  (mils) 

5.  The  program  METRIC  may  be  written  on  a  tape  which  is  used  in 
each  run  (see  section  on  control  cards). 

Intermediate  tapes  Note:  On  the  UN  I  VAC  110R,  intermediate  data  may 
he  written  on  magnetic  drum  rather  than  magnetic  tape. 

1.  When  a  UNI!.0ti  tape  is  being  processed  and  it  contains  data 
from  two  radars,  the  data  from  both  radars  may  be  reduced  in  a  single 
run.  In  this  case,  when  the  data  from  the  first  radar  is  being  pro¬ 
cessed,  information  from  the  second  radar  is  being  written  on  Fortran 
unit  i)  in  an  unformatted  mode. 

Word 

1  Time  (seconds  from  1 i f t-of f - t ime ) 

2  Slant  range  from  radar  (meters) 

ii  Azimuth  from  radar  (radians) 

A  Elevation  from  radar  (radians) 
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2  .  In  any  run  with  wind  data  to  p;.  ;  I'ssni,  the  miiivit1ii.il 

wind  samples  art'  ivt  orded  on  Korlr.it-.  i^rital  unit  1  5  tor  use  in 
obtaining  averaged  layer  wind  data.  It  the  wind  data  is  input  >ui 
(aiats,  the  input  itata  points  are  returded.  When  tlie  input  is  digital 
tape,  data  is  computed  and  ston'd  at  one  sample  per  scmnd.  Iliese 
rercirtis  are  written  in  an  unformatted  mode. 


Wo  rd 

1  Time  (seconds  from  lift-off  time) 

2  Altitude  (meters) 

.1  Alt  it  ude  (  feet ) 

4  N-S  wind  component  (meters/sec) 

5  K-V  wind  component  (meters/sec) 

t>  Wind  speed  (meters/sec) 

7  Wind  direction  (degrees) 

8  Kail  velocity  (meters/sec) 

11  Number  of  points  since  last  good  point  which  were  rejected 
as  had  data. 

Output  tapes  At  present  there  art'  no  output  tapes  being  gecierat< 
by  MKTK0C . 


OUTPUT  LISTINGS 

The  output  listings  are  in  metric  units,  with  additional  listings 
in  English  units  available  as  an  option.  The  data  consists  ot  : 

1.  Wind  Data 

A.  Individual  wind  samples 

1.  One  sample/second  with  digital  tape  input 

2.  Input  points  with  punch  card  input 
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B.  Averaged  wind  data  for  even  kilometer  levels  average d 

over  2  kilometer  layers. (When  temperatures  are  included, 
the  averaged  winds  are  listed  combined  with  the  thermo¬ 
dynamic  data.) 

II.  Thermodynamic  Data 

A.  Type  of  data 

1.  Wind  averaged  over  2  kilometer  layers 

2.  Temperatures 

3.  Pressure  computed  using  hydrostatic  ctpial ions , 
rocket  temperatures,  and  base-level  point  1  nun 
conjunctive  r.iob 

4.  Computed  density 

5.  Computed  speed  of  sound 

B.  Types  of  levels  for  thermodynamic  data 

1.  Significant  temperature  levels 

2.  liven  kilometer  or  .5  kilometer  levels 

3.  5000  fool  or  1000  foot  levels 

4.  Mandatory  (constant-pressure)  levels 

C.  Types  of  listings  output 

1 .  Standard  data  summary 

2.  MKN  formatted  data 

3.  Ozonesonde  parameters  (optional) 

4.  English  units  data  (optional) 
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XII  FLOWCHARTS  OF  OVERALL  SYSTEM  OPERATION 


i 


i 


i 


r 


Initialization 
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nr!  Data  Card  Processing 


Print 

heading 


Compute  km  level 
altitudes, 
interpolate  for 
tempo ratu res 


.T.essa.^e 


XL  I  1  FLOWCHARTS  OF  MODIFIED  OR  NEW  PROGRAMS 


FLOWCHART  FOR  SUBIIO  L'T  I N  E  RT  DAT  A  -  2.1 
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FLOWCHART  FOR  SUBROUTINE  RTDATA  -  1. 
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